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Divergences of the localization lengths in the two-dimensional, off-diagonal Anderson 

model on bipartite lattices 
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1. Introduction 

The disorder-induced metal-insulator transition (MIT) 
and the concept of Anderson localization^"^) has been the 
subject of intense research activities for more than forty 
years. The scaling theory of localization^^ predicts that 
for generic, interaction-free situations in 2D all states 
are localized; there is no disorder-driven MIT^^^^ and 
the system remains an insulator. However, non-localized 
states were found at the band center of an Anderson 
model with purely off-diagonal disorder. ^""^'^^ Numer- 
cal studies of such systems revealed that the localization 
length diverges at the energy E = q}5,16} Scaling prop- 
erties of this divergence suggest that the states at the 
band center are critical^®' and that are indeed no truly 
extended states in 2D even in this off-diagonal case. 

Of special importance is the model on a bipartite lat- 
tice with an even number of sites where the energy spec- 
trum is strictly symmetric around E = 0. In this case, 
related to the chiral universality class, E — states are 
non-localized in any dimension. ^^"^"^ This has been re- 
cently demonstrated analytically in 2D and 3D using a 
renormalization group (RG) approach.2i) In Ref.^i) it 
has been also suggested that for sufficiently large ener- 
gies the divergence of the localization length at the band 
center may be described by a power law. 
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whereas it takes more complicated form 
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below a certain crossover energy E* . 

Recently^^) we studied such 2D bipartite systems with 
various types of off-diagonal disorder by means of the 
transfer-matrix method (TMM) and investigated the di- 
vergence of the localization length close to E — 0. We 
showed that the data fit the power-low behaviour down 
to energy E — 2 x 10~^ and the corresponding critical 
exponents seem to depend on the type of disorder dis- 
tribution and strength. Here we extend these results to 
even smaller energies and examine honeycomb lattices 
where the Van Hove singularity at E — does not inter- 
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fere with the divergence due to the bipartiteness. 

2. Computation of the localization lenghts and 
the critical exponents 

We considered a single electron on a 2D lattce de- 
scribed by the Hamiltonian 
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where \i) denotes the electron at site i. The off-diagonal 
disorder was introduced by random hopping elements tij 
between nearest neighbor sites. We considered bipartite 
square and honeycomb lattices. The latter is topologi- 
cally equivalent to a brick-layer structure, ^'^^ which im- 
plies that one of the connections to the nearest-neighbor 
is absent — the corresponding hopping element is equal 
to 0. 

On both lattices we examined three different distribu- 
tions of the off-diagonal elements tij : 
(1) a rectangular distribution^^) 
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if \Uj - c| < w/2, 
otherwise, 



(2) a Gaussian distribution 
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(3) a rectangular distribution of the logarithm of tij^"^^ 



1/w if llnty/iol < W2, 
otherwise. 



In the case of rectangular and Gaussian distributions of 
hopping elements we have set the width w and the stan- 
dard deviation a of the distribution to 1 and centered 
it at c = 0. In the case of logarithmic-t distribution we 
chose = 1 which set the energy scale and performed 
the calculations for two values of the distribution width, 
e.g., w ^2 and 6. 

We used the transfer-matrix method^'^'^s) to compute 
the localization lengths for strips of widths M up to 220 
in the energy interval 10~^ < E < 0.1, the actual values 
depend on the disorder parameters. The accuracy of the 
results was 1%. Next, finite-size scaling^^) was applied 
to the data. We used M as large as possible in order 



2 



Short Note 



to avoid problems related to the finite-size effects. For 
some cases this imposes the restriction on the smallest 
energies which we were able to use in FSS. 

To investigate the divergence of the localization 
lenghts at the band center we plotted the scaling pa- 
rameter ^ in a double-logarithmic plot. The deviation 
of the divergence from the power-law behaviour should 
then be easily seen. In most cases we observed that at 
the energies close to E = the divergence is slower than 
described by a power-law. However, these deviations get 
smaller when we increase widths M of the systems used 
in FSS (an example of such behaviour is shown in Fig. 
1), therefore we relate them to finite-size effects. 
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Fig. 1. Scaling parameter g vs. energy. Left panel: square lattice, 
Gaussian distribution (ct = 1, c = 0), right panel: honeycomb 
lattice, logarithmic t distribution (ui = 6); filled symbols: results 
for AI = 50-100, open symbols: results for M = 110-160 (Gaus- 
sian distribution) or M = 100-150 (logarithmic distribution). 



3. Results and disussions 

In Table I we collected the critical exponents obtained 
for different lattices and disorder distributions. All val- 
ues are in the range of 0.2 — 0.6 and depend on the disor- 
der type. Generally the exponents are larger for weaker 
disorders (i.e. when the localization lengths are larger). 
The differences between square and honeycomb lattice in 
most cases do not exceed the accuracy of the results. The 
exception is the logarithmic t distribution with w = 2, 
when the exponent is almost 2 times larger on a hon- 
eycomb lattice. The large change in the exponent may 
be related to the low density of states on a honeycomb 
lattice, where the van Hove singularity in the density of 
states is absent. 

The localization lengths exhibit power-law behaviour 
in a wide energy range, the lower lower bound ii'min of 
which is indicated in Table I. For smaller energies we ob- 
serve some deviations, however, these become smaller for 
larger system widths M we would argue that these are 
finite-size effects that will eventually vanish if one uses 
large enough system sizes. The exception to this be- 
haviour is the logarithmic-i disribution for w = 6, where 
below an energy E !v 10~^ the deviation from power law 
seems to be size independent. This may be an indica- 
tion of the crossover to the non-power-law behavior as 
predicted in Ref.^^-' or to a power-law with a different 



Table I. Estimated values of the exponents u of the localization 
lengths for various disorder strengths and distributions. The 
range of the strip widths M is displayed in the second column; 
the power-law behaviour is obeyed in the energy ranging from 
up to 1 X 10 ^. The error bars represent the standard 
deviations from the power-law fit and should be increased up to 
one order of magnitude for a reliable representation of the actual 
errors. 

square lattice 



disorder 
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^min 
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box, c = 




150 - 220 


8 X 10" 


7 


0.317 ±0.007 


Gaussian, 


c = 


110 - 160 


2 X 10" 


6 


0.303 ± 0.006 


Int, w = 2 




110 - 190 


2 X 10- 


6 


0.357 ±0.009 


Int, w = 6 




120 - 170 


1 X 10" 


5 


0.232 ± 0.007 


honeycomb lattice 


disorder 




M 


-^min 




V 


box, c = 




110 - 170 


2 X 10" 


8 


0.290 ±0.004 


Gaussian, 


c = 


120 - 170 


1 X 10" 


6 


0.273 ±0.005 


Int, w = 2 




100 - 170 


2 X 10" 


4 


0.604 ± 0.015 


Int, w = a 




120 - 170 


1 X 10" 


5 


0.238 ± 0.007 



exponent. However, there is also a possibilty that this 
may be an effect of numerical problems which appear for 
strong logarithmic disorders. 
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